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INTRODUCTION 



This subroutine set provides an IMP-16 microprocessor with double-precision and floating-point capability. 
Double-precision and floating-point computation is very useful when additional precision is needed or the 
range of numbers is expected to be large. In the following pages, individual subroutines are described along 
with necessary user instructions. Subroutines included in the set are listed in the table below. 
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GENERAL CHARACTERISTICS 



The arithmetic subroutine set uses 512 memory words located at addresses FCOO through FDFF (64512 
through 65023, decimal). All subroutines are contained wholly within the set; they require no external 
program code. The IMP-16 extended instruction set is required. Operands and results, in most cases, 
use the general registers. The only memory locations written into by the set are in the range 00E0 
through 00EF. All constants are contained within the subroutine set. 

EXECUTION TIMES 



Execution times shown on each of the subroutine data sheets are expressed in terms of three variables, where 

R = Number of main memory read cycles 
W = Number of main memory write cycles 
N = Number of microprogram cycles 
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These execution time expressions can thus be used to calculate subroutine execution time for any IMP-16 
implementation. A few subroutine execution times are also dependent on user -specified variables (such as 
shift counts in shift instructions). These are defined on the individual data sheets. All execution times are 
based on worst -case conditions. 



NUMBER REPRESENTATION 



Numbers in a computer program often represent real world parameters. In this subroutine set, four different 
modes of numerical representation are used: integer, fractional, double-precision, and floating-point. All 
four modes have important applications in IMP-16 systems. 

Integer notation is the most common and obvious mode of numerical representation. It is simply an integer 
count of something. In the IMP-16, numbers up to 2 16 may be represented because of the 16-bit word length. 
To represent both positive and negative integers, a twos -complement system is used. This gives an integer 
range of 

-2 15 < integer < 2 15 

With twos -complement notation, the most significant bit, bit 15, can be thought of as a sign bit, and the binary 
point can be thought of as being to the right of bit 0. 



FRACTIONAL NOTATION 



Many real -world parameters are continuous rather than discrete. Such values can never be indicated with per- 
fect precision because of the limited word length of any computer. Take the airspeed of an aircraft as an 
example. It may have a potential range from zero to 1000 miles per hour. We would rarely expect its speed 

to be an exact round number. At some instant the speed may be 556. 234 If we want to express this 

speed as precisely as possible given a 16-bit word length, an integer notation of 556 miles per hour would not 
be the best choice. A more convenient technique is to use fractional notation. The parameter, in this case 

speed, can be scaled so that 1000 miles per hour is equivalent to 1. 0000 This would allow the full 16-bit 

capacity to be utilized. After scaling, the range of numerical expression of a fractional number is 

- 1.0 < N < 1.0 

With fractional notation, the binary point can be thought of as being to the right of bit 15. Thus, the largest 
positive binary number that can be represented is 0. 11111. . . or almost one. Precision is about 3.1* 10~ 5 . 



DOUBLE PRECISION 



When precision greater than 16-bits is required, double -precis ion numbers can be used. Double-precision 
numbers give a precision of one part in 2 31 or about 4. 6* 10" 10 . Double-precision numbers are an extended 
fractional notation. Note that two IMP-16 words are required. 



FLOATING POINT 



Floating-point representation is a technique to express numbers in a form similar to scientific notation, with 
a fractional mantissa and an exponent. The value of any floating-point number so represented is equivalent to 
mantissa * (2**exponent). The mantissa consists of 24 bits and gives a precision of one part in 2^ 3 . The 
exponent has 8 bits and allows exponents in the range 

7 7 
-2 < exponent < 2 
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The great advantage of floating-point numbers is that the scaling factor of a fractional number becomes a 
visible attribute which can be manipulated. Thus, the scaling factor can be altered to maintain precision 
or facilitate computations. 



NUMBER REPRESENTATION EXAMPLES 
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• INTEGER 



15 



-tt- 



0 



BINARY POINT 



EXAMPLES: 

5 = 0005i6 
-5 = FFFB 16 



RANGE: -32768 < N < 32768 
MAXIMUM PRECISION: 1 
SMALLEST VALUE: 1 



SPEED = 556.234 = 556 =0220^ 



• FRACTIONAL 



15 . 14 



\ 



BINARY POINT 
RANGE: -1 < N < 1 

MAXIMUM PRECISION: l/2 15 = 3.1 * 10 " 5 
SMALLEST VALUE: 2" 15 



EXAMPLES: 
1/2 = 4000 16 
-1/2 =C000 16 
-1 = 8000 16 

556.234.. 
SPEED = 1000.000 



= 0. 556234 10 = 4732 



s 16 



DOUBLE PRECISION (FRACTIONAL) 



31. 30 



RANGE: -1< N < 1 

MAXIMUM PRECISION: l/2 31 = 4. 6 * 10~ 10 
SMALLEST VALUE: 2~ 31 
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15 



EXAMPLES: 
1/2 = 4000 0000 16 
-1/2 =C000 0000 16 
-1 = 8000 0000, 



'16 



SPEED = 4732 ACFB^ 



• FLOATING POINT 



RANGE: -2 127 < N< 2 127 
PRECISION: 2~ 23 ( 2 EXPONENT ) 
SMALLEST MAGNITUDE: 2~ 151 
LARGEST MAGNITUDE: (l-2~ 23 ) 2 127 



23. 22 





1 — 11 — 1 


1 n 1 






7 ».. 0 



MANTISSA (M) 



N = M * 2 



E 



EXPONENT (E) 



EXAMPLES: 
1/2 = 4000 00i 6 (00i6) 
1/2 = 2000 OOie (01i 6 ) 
1/2 = 1000 00i6 (02i 6 ) 
1/4 = 2000 00i 6 (00i 6 ) 
1/4 = 4000 00i6 (FFi 6 ) 



SPEED = 556. 234 = 4587 7A i6 (0A 16 ) 
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ANGULAR REPRESENTATION 

Angles are expressed as a double -precision fraction of tr radians. 




NOTATION 

In the descriptions which follow, RO, Rl, R2, and R3 represent the contents of the four accumulators = AGO, AC1, 
AC2, and AC3. 

Normalization of a floating point number means that the mantissa is shifted left until the most significant bit is a 1, 
and the exponent adjusted accordingly. 
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SUBROUTINE: SINGLE PRECISION MULTIPLY 



LABEL: MULT 



DESCRIPTION: MULTIPLIES TWO SINGLE -PRECISION, SIGNED, FRACTIONAL 16-BIT NUMBERS TO 
GIVE A SIGN PLUS A 30-BIT RESULT 



ENTRY POINT: FC01 

ENTRY CONDITIONS: R Q = X R 2 = Y 

EXIT CONDITIONS: Rq, R^ = X * Y R 2 - ALTERED R g = UNCHANGED 

EXECUTION TIME: 30R + 243N + W 



TEMPORARY LOCATIONS USED: 00E0 

SUBROUTINES CALLED: DSHL, DPCOMP 

STACK WORDS NEEDED: 2 
(MAXIMUM DEPTH) 



COMMENTS: 1. ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

2. FRACTIONAL DIFFERS FROM AN INTEGER MULTIPLY IN THAT THE RESULT OF A 
FRACTIONAL MULTDPLY IS SHIFTED LEFT 1 BIT. 



SUBROUTINE: SINGLE PRECISION DIVIDE 



LABEL: DW 



DESCRIPTION: DIVIDES A SIGNED, FRACTIONAL 16-BIT NUMBER INTO A SIGNED, FRACTIONAL 
32-BIT NUMBER, PROVIDING A QUOTIENT AND AN INTEGER REMAINDER 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 
EXECUTION TIME: 



FC1A 



Rq, R-l = DIVIDEND 
R 0 = QUOTIENT 
49R + W + 343N 



R 2 = DP7ISOR 



R x = REMAINDER R 2 = ALTERED R 3 = UNCHANGED 



TEMPORARY LOCATIONS USED: 00E0 

SUBROUTINES CALLED: DSHR, DPCOMP 

STACK WORDS NEEDED: 2 
(MAXIMUM DEPTH) 



COMMENTS: 1. ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

2. OVERFLOW FLAG IS SET IF MSB OF DP7IDEND > DIVISOR. 
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SUBROUTINE: DOUBLE PRECISION MULTIPLY 



LABEL: DPMUL 



DESCRIPTION: MULTIPLIES DOUBLE -PRECISION OPERANDS TO GIVE A DOUBLE -PRECISION PRO- 
DUCT 



ENTRY POINT: FCCF 

ENTRY CONDITIONS: R Q , R x = X R 2 , R 3 = Y 

EXIT CONDITIONS: R Q , R-^ = X * Y R 2 , R 3 ALTERED 

EXECUTION TIME: 200R + 8W + 1125N 



TEMPORARY LOCATIONS USED: 
SUBROUTINES CALLED: 
STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



00E0, 00E1, OOEE, OOEF 
QUAD, MULT, DPCOMP, DSHL 
4 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: POSITIVE DOUBLE PRECISION DIVIDE 



LABEL: DPDF7 



DESCRIPTION: UNSIGNED FRACTIONAL DIVIDE OF R 2 , R 3 INTO R Q , R p BOTH DOUBLE -PRECISION 
POSHTVE FRACTIONS 



ENTRY POINT: FDE3 

ENTRY CONDITIONS: R 0 , R 1 = POSITIVE DIVIDEND R 2 , R g = POSITIVE DIVISOR 

EXIT CONDITIONS: R 0 , R x = QUOTIENT R 2 , R 3 = REMAINDER SELFF = 0 

EXECUTION TIME: 510R + 97W = 2578N 



TEMPORARY LOCATIONS USED: 00E0, 00EC-00EF 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 0 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: DOUBLE PRECISION SQUARE 



LABEL: DPSQUARE 



DESCRIPTION: SQUARES THE DOUBLE -PRECISION NUMBER IN R Q , R-^ 



ENTRY POINT: FCCD 

ENTRY CONDITIONS: R Q , R x = X 

EXIT CONDITIONS: R Q , Rj = X 2 R 2 , R3 ALTERED 

EXECUTION TIME: 202R + 8W + 1137N 



TEMPORARY LOCATIONS USED: 00E0, 00E1, OOEE, OOEF 

SUBROUTINES CALLED: DPMULT, QUAD, MULT, DPCOMP, DSHL 

STACK WORDS NEEDED: 4 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSED ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: DOUBLE PRECISION COMPLEMENT 



LABEL: DPCOMP 



DESCRIPTION: COMPUTES TWOS COMPLEMENT OF A DOUBLE -PRECISION NUMBER 



ENTRY POINT: FDBA 

ENTRY CONDITIONS: R Q , Rj = X 

EXIT CONDITIONS: R 0 , Rj = TWOS COMPLEMENT OF X R 2 , R3 PRESERVED SELFF = 0 

EXECUTION TIME: 16R + 53N 



TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 0 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: DOUBLE PRECISION ARITHMETIC SHIFT 



LABEL: DPSH 



DESCRIPTION: SHIFTS A DOUBLE -PRECISION NUMBER TO RIGHT OR LEFT AS SPECIFIED BY R 2 . 



ABSOLUTE VALUE OF R 2 INDICATES NUMBER OF BITS TO SHIFT. 
INDICATES LEFT SHIFT; NEGATIVE R 2 MEANS RIGHT SHIFT. 



POSITIVE R. 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 
EXECUTION TIME: 



FDDB 

R 2 = M (SHIFT COUNT) 
SHIFTED X R 2 = ALTERED 
LEFT X 5R + 20N + [8R + 32N] M 
RIGHT 7R + 26N + [4R + 21N] M 



Rq> Rj_ - X 
Rq, Rj 



R 3 PRESERVED 



TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 0 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

RIGHT SHIFT IS ARITHMETIC; THAT IS, SIGN BIT IS PRESERVED 



SUBROUTINE: DOUBLE PRECISION ARITHMETIC SHIFT RIGHT 



LABEL: DPSHR 



DESCRIPTION: SHIFTS DOUBLE-PRECISION NUMBER RIGHT M BITS, WHERE R 2 = M 



ENTRY POINT: FDCE 

ENTRY CONDITIONS: Rq, Rj = X R 2 = M (BITS TO SHIFT) 

EXIT CONDITIONS: R 0 , R x = RIGHT SHIFTED X R 2 ALTERED Rg PRESERVED 

EXECUTION TIME: 5R + 20N + (4R + 21N) M 



TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 0 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
SIGN BIT IS PRESERVED 
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SUBROUTINE: DOUBLE PRECISION SHIFT LEFT 



LABEL: DPSHL 



DESCRIPTION: SHIFTS DOUBLE -PRECISION NUMBER LEFT M BITS, WHERE R 2 = M 



ENTRY POINT: FDDB 

ENTRY CONDITIONS: R Q , R 1 =X R g = M (NUMBER OF BITS TO SHIFT) 

EXIT CONDITIONS: Rq, R ± = LEFT SHIFTED X R 2 ALTERED R 3 PRESERVED 

EXECUTION TIME: 5R + 20N + (8R + 32N) M 



TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 0 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: QUADRANT TESTS 



LABEL: QUAD 



DESCRIPTION: GIVEN TWO SIGNED DOUBLE -PRECISION NUMBERS CI AND C2, THE INITIAL SIGN OF 
EACH IS INDICATED BY A FLAG IN MEMORY AND THE ABSOLUTE VALUE OF THE 
NUMBERS IS GIVEN. 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 



EXECUTION TIME: 



FDAO 
R 



R l ~°1 



0' "1 1 —'£> 
R Q , R-l = ABS (C x ) 



47R + W + 175N 



R 3 ~ C 2 



R 2 , R 3 = ABS (C 2 ) 



00E1 = 0 + 
1 + 
2 
3 

WHERE ZERO IS + 



TEMPORARY LOCATIONS USED: 00E1 

SUBROUTINES CALLED: DPCOMP 

STACK WORDS NEEDED: 1 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: SINE 



LABEL: SIN 



DESCRIPTION: CALCULATES SINE OF DOUBLE -PRECISION ANGLE IN R 0 , Rj 
CALCULATION IS MADE: SIN(X) = COS (X - Or/ 2) 
ACCURACY IS 0. 00000 0002^ 1Q j 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 
EXECUTION TIME: 



FCFA 

Rq> r^ = x 

R 0 , R x = SIN(X) R 2 , R 3 ALTERED 
1697R + 82W + 9300 N 



TEMPORARY LOCATIONS USED: 
SUBROUTINES CALLED: 
STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



00E0-00EB, 00EE, 00EF 
DPCOMP, DPMULT, PEXPN 
6 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: COSINE 



LABEL: COS 



DESCRIPTION: CALCULATES COSINE OF DOUBLE -PRECISION ANGLE IN R Q , Rj 
POLYNOMIAL APPROXIMATION USED IS: 

COS(X) = 1 + X 2 (A 2 + X 2 (A 4 + X 2 (A 6 + X 2 (A g + X 2 (A 1Q ))))) 
WHERE: A 2 = -0.49999 99963 A g = 0. 00002 47609 

A 4 = 0. 04166 66418 A 10 = -0. 00000 02605 

A 6 = -0. 00138 88397 

FOR ANGLES IN THE FIRST QUADRANT 
ACCURACY IS 0. 00000 0002^ 10) 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 
EXECUTION TIME: 



FCFC 

R 0 , Rj = COS(X) R 2 , Rg ALTERED 
1677R + 82W + 9285N 



TEMPORARY LOCATIONS USED: 0OEO-00EB, 00EE, 00EF 

SUBROUTINES CALLED: DPCOMP, DPMULT, PEXPN 

STACK WORDS NEEDED: 6 
(MAXIMUM DEPTH) 



COMMENTS: 



ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: ARCTANGENT 



LABEL: ARC TAN 



DESCRIPTION: COMPUTES ARCTANGENT (C1/C2), WHERE CI AND C2 ARE DOUBLE -PRECISION 
NUMBERS IN THE RANGE: 

-1 < C1,C2 < 1 

THE POLYNOMIAL APPROXIMATION USED WAS: 



ARCTAN (X) = X (1 + X 2 (A 2 + X 2 (A 4 + X 2 (A 6 + X 2 (A 8 + X 2 (A 10 + X 2 (A 12 + X 

(A 14 + X 2 (A 16 )))))))) WHERE 

A 2 = -0. 33333 14528 A 10 = -0. 07528 96400 

A 4 = 0.19993 55085 A 12 = 0-04290 96138 

A 6 = -0.14208 89944 A 14 = -0. 01616 57367 

A g = 0. 10656 26393 A 16 = 0. 00286 62257 

ACCURACY IS 0.00000 002(10) 
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ENTRY POINT: FD37 

ENTRY CONDITIONS: R 0 , Rj = C 1 R 2 , Rg = C 2 

EXIT CONDITIONS: R Q , R 1 = ARCTAN (C-j/C^ Ro, Rg ALTERED 

EXECUTION TIME: 2985R + 231W + 15892N 



TEMPORARY LOCATIONS USED: 00E0-00EF 

SUBROUTINES CALLED: QUAD, DPCOMP, DPDF7, DPMUL, PEXPN, DSHR 

STACK WORDS NEEDED: 7 
(MAXIMUM DEPTH) 



COMMENTS: 1. ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

2. SINGULARITY CASE WHERE CI = C2 = 0 YIELDS ARCTAN(Cl/C2) = 0. 



SUBROUTINE: FLOATING POINT ADD 



LABEL: FPADD 



DESCRIPTION: ADDS TWO FLOATING-POINT NUMBERS. 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 
EXECUTION TIME: 



FC30 
Rnj Hi 



= X 



R, 



R Q> R x = X + Y 
147R + 7W + 613N 



2' 



R 3 =Y 



l 2> 



R 3 ALTERED 



TEMPORARY LOCATIONS USED: 
SUBROUTINES CALLED: 
STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



00E0, 00E1, 00E2 

EXTEXP, DSHR, CZERO, ADD EXP 
1 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: FLOATING POINT MULTIPLY 



LABEL: FPMUL 



DESCRIPTION: MULTIPLIES TWO FLOATING-POINT NUMBERS 



ENTRY POINT: FC5D 

ENTRY CONDITIONS: R Q , n 1 =X R 2> R3 = Y 

EXIT CONDITIONS: R Q , R^ = X * Y R 2 , Rg ALTERED 

EXECUTION TIME: 1215R + 63W + 5077N 



TEMPORARY LOCATIONS USED: 00E0, 00E1, 00E5, 00E6, OOEE, OOEF 

SUBROUTINES CALLED: DLNORM, DPMUL, ADD EXP, CZERO 

STACK WORDS NEEDED: 4 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: POSITIVE FLOATING POINT DIVIDE 



LABEL: FPDIV 



DESCRIPTION: DP7IDES TWO POSITIVE FLOATING-POINT NUMBERS 



ENTRY POINT: FC6C 

ENTRY CONDITIONS: R Q , R ± = X R 2 , Rg = Y 

EXIT CONDITIONS: R Q , R 1 = X/Y 

EXECUTION TIME: 1540R + 152W + 6584N 



TEMPORARY LOCATIONS USED: 00E0, 00E1, 00E2, 00E5, 00E6, OOEC, OOED, OOEE, OOEF 

SUBROUTINES CALLED: DLNORM, DSHR, DPDP7, ADDEXP, CZERO 

STACK WORDS NEEDED: 4 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: FLOATING POINT COMPLEMENT 



LABEL: FPCOMP 



DESCRIPTION: COMPLEMENTS A FLOATING-POINT NUMBER IN R Q , R ± 



ENTRY POINT: FCB9 

ENTRY CONDITIONS: R Q = X x R 1 = X 2 , EX 

EXIT CONDITIONS: R Q = Y{ R-, = X^, EX R 2 , R 3 PRESERVED 

EXECUTION TIME: 50R + 2W + 183N 



TEMPORARY LOCATIONS USED: 00E0 

SUBROUTINES CALLED: EXTEXP, DPCOMP, ADDEXP 

STACK WORDS NEEDED: 1 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: CHECK ZERO EXPONENT 



LABEL: CZERO 



DESCRIPTION: IF MANTISSA OF A FLOATING-POINT NUMBER IS ZERO, EXPONENT IS FORCED TO 
ZERO. 



ENTRY POINT: FC63 
ENTRY CONDITIONS: R 0 , R ± = X ls X 2 (EX) 



EXIT CONDITIONS: IF X ls X 2 = 0 R Q , Rj_ = 0,0(0) ^ R PRESERVED 

IF X lf X 2 f 0 R Q , R 2 = X v X 2 (EX) J 



EXECUTION TIME: 36R + 2W + 139N 



TEMPORARY LOCATIONS USED: 00E0 

SUBROUTINES CALLED: EXTEXP, ADDEXP 

STACK WORDS NEEDED: 1 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: EXTRACT EXPONENT TO STACK 



LABEL: EXTEXP 



DESCRIPTION: EXTRACTS EXPONENT FROM A FLOATING-POINT NUMBER, STORES EXPONENT ON 
STACK, AND CLEARS EXPONENT FIELD OF NUMBER. 



ENTRY POINT: FC7C 

ENTRY CONDITIONS: R Q , R ± = X v X 2 (EX) 

EXIT CONDITIONS: R Q , R x = X 1§ X 2 (0) STACK = EX 

EXECUTION TIME: 13R + W + 51N 



R 2 , R g PRESERVED 



TEMPORARY LOCATIONS USED: 00E0 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 0 
(MAXIMUM DEPTH) 



COMMENTS: 



ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: ADD EXPONENT FROM STACK 



LABEL: ADDEXP 



DESCRIPTION: AN EXPONENT IN THE STACK IS PLACED IN A FLOATING-POINT WORD IN R Q , R r 



R 0 



Ri 



STACK 







x 2 


x 3 


BEFORE 








X l 




X 2 


EX 




AFTER 



ENTRY POINT: FCBD 

ENTRY CONDITIONS: R Q , ^ = X p X 2 (X 3 ) STACK = EX 
EXIT CONDITIONS: R Q , R 1 = X r X 2 (EX) R 2 , Rg PRESERVED 



EXECUTION TIME: 



13R + W + 51N 



TEMPORARY LOCATIONS USED: N 00E0 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: 1 FOR INPUT EXPONENT ONLY 
(MAXIMUM DEPTH) 



COMMENTS: 1. ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

2. WARNING: ACCESS THIS SUBROUTINE ONLY BY A JSR INSTRUCTION! 
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SUBROUTINE: LEFT NORMALIZE 



LABEL: LFNOR 



DESCRIPTION: LEFT NORMALIZES A FLOATING-POINT NUMBER IN R Q , Rj 



ENTRY POINT: 
ENTRY CONDITIONS: 
EXIT CONDITIONS: 
EXECUTION TIME: 



FC97 



R 



0' 



Rj =X 



R 0 , R x = X NORMALIZED 
449R + 24W + 1745N 



;R 2 , R3 ARE PRESERVED 



TEMPORARY LOCATIONS USED: 
SUBROUTINES CALLED: 
STACK WORDS NEEDED: 
(MAXIMUM DEPTH) 



00E0, 00E1 

CZERO, FPCOMP, EXTEXP, DSHL, ADDEXP 
2 



COMMENTS: 



ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: DOUBLE LEFT NORMALIZE 



LABEL: DLNORM 



DESCRIPTION: NORMALIZES TWO FLOATING-POINT NUMBERS R Q , R x AND R 2 , Rg, WITH THE 

RESULTANTS SWITCHED (IN REGISTERS R 2 , Rg AND R Q , R x , RESPECTIVELY). NEW 
EXPONENTS ARE IN MEMORY. 



ENTRY POINT: FC8A 

ENTRY CONDITIONS: R Q , R x = Xj, X 2 (EX) Rg, Rg = Y v Y £ (EY) 

EXIT CONDITIONS: R 0 , R^ = Y NORMALIZED (EX = 0) R 2 , Rg = X NORMALIZED (EX = 0) 

00E6 = EY 00E5 = EX 

EXECUTION TIME: 947R + 52W + 3691N 



TEMPORARY LOCATIONS USED: 00E0-00E2, 00E5, 00E6 

SUBROUTINES CALLED: LFNOR, EXTEXP 

STACK WORDS NEEDED: 3 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 
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SUBROUTINE: FRACTIONAL TO FLOATING POINT 



LABEL: FLOAT 



DESCRIPTION: CONVERTS DOUBLE -PRECISION FRACTIONAL NUMBER TO FLOATING-POINT 
FORMAT BY TRUNCATING LEAST SIGNIFICANT 8 BITS. 



ENTRY POINT: 


FC85 








ENTRY CONDITIONS: 


Rq> r^ = 








EXIT CONDITIONS: 


R 0 , R-,^ = 


Xi, 


X 2 « (0) 


R 2 , R 3 PRESERVED 


EXECUTION TIME: 


3R + 9N 







TEMPORARY LOCATIONS USED: NONE 

SUBROUTINES CALLED: NONE 

STACK WORDS NEEDED: NONE 
(MAXIMUM DEPTH) 



COMMENTS: 



ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 



SUBROUTINE: FLOATING POINT TO FRACTIONAL 



LABEL: SFO 



DESCRIPTION: CONVERTS A FLOATING-POINT NUMBER TO A FRACTIONAL DOUBLE -PRECISION 
NUMBER 



R 0 



Ri 



EX 



BEFORE 



MANTISSA 



EXPONENT 



X 



AFTER 



DOUBLE PRECISION 



ENTRY POINT: 


FCC7 




ENTRY CONDITIONS: 


R 0 , R x = X X> X 2 (EX) 




EXIT CONDITIONS: 


R 0 , R 1 = X DOUBLE PREC. 


R 2 , R 3 PRESERVED 


EXECUTION TIME: 


207R + W + 822N 





TEMPORARY LOCATIONS USED: 00E0 

SUBROUTINES CALLED: EXTEXP, DSH 

STACK WORDS NEEDED: 2 
(MAXIMUM DEPTH) 



COMMENTS: ALL ADDRESSES ARE EXPRESSED IN HEXADECIMAL NOTATION. 

CAUTION: VALUE OF INPUT MUST BE < 1. 0 IN MAGNITUDE. 
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National Semiconductor Corporation 

2900 Semiconductor Drive 
Santa Clara, California 95051 
(408) 732-5000 
TWX: 910-339-9240 



National Semiconductor GmbH 

D 808 Fuerstenfeldbruck 
Industriestrasse 10 
West Germany 
Telephone: (08141) 1371 
Telex: 05-27649 



NS Electronics SDN BHD 

Batu Berendam 
Free Trade Zone 
Malacca, Malaysia 
Telephone: 5171 

Telex: NSELECT 519 MALACCA (c/o Kuala Lumpur) 



National Semiconductor (UK) Ltd. 

Larkfield Industrial Estate 
Greenock, Scotland 
Telephone: GOUROCK 33251 
Telex: 778 632 



NS Electronics (PTE) Ltd. 

No. 1100 Lower Delta Rd. 
Singapore 3 
Telephone: 630011 
Telex: NATSEMI RS 21402 



REGIONAL AND DISTRICT SALES OFFICES 



ALABAMA 

DIXIE REGIONAL OFFICE 
3322 Memorial Parkway, S.W. #67 
Hunlsville, Alabama 35802 
(205) 881-0622 
TWX: 810-726-2207 

ARIZONA 

♦ROCKY MOUNTAIN REGIONAL OFFICE 
7349 Sixth Avenue 
Scotlsdale, Arizona 85251 
(602) 945-8473 
TWX: 910-950-1195 

CALIFORNIA 

* NORTH-WEST REGIONAL OFFICE 
2680 Bayshore Frontage Road, Suite 112 
/Mountain View, California 94043 

(415) 961-4740 
TWX: 910-379-6432 

NATIONAL.SEMICONDUCTOR 
*LOS ANGELES REGIONAL OFFICE 

Valley Freeway Center Building 

15300 Ventura Boulevard, Suite 305 

Sherman Oaks, California 91403 

(213) 783-8272 

TWX: 910-495-1773 

NATIONAL SEMICONDUCTOR 
♦SOUTHERN CALIFORNIA REGIONAL OFFICE 

17452 Irvine Boulevard, Suite M 

Tustin, California 92680 

(714) 832-8113 

TWX: 910-595-1523 

CONNECTICUT 

DISTRICT OFFICE 
Commerce Park 
Danbury, Connecticut 06810 
(203) 744-2350 
♦MID-ATLANTIC REGIONAL SALES OFFICE 
25 Sylvan Road South 
Westport, Connecticut 06880 
(203) 226-6833 

INTERNATIONAL SALES OFFICES 
AUSTRALIA 

NS ELECTRONICS PTY, LTD. 
Cnr. Stud Road & Mountain Highway 
Bayswater, Victoria 3153 
Australia 

Telephone: 03-729-6333 
Telex: 32096 

BELGIUM 

NATIONAL SEMICONDUCTOR BELGIUM 
789 Ave. Houba de Strooper 
1020 Bruxelles 
Telephone: 02-478-3400 
Telex: 61 007 NatSem B 

CANADA 

NATIONAL SEMICONDUCTOR 
DISTRICT OFFICE 
1111 Finch Avenue West 
Suite 154 

Downview, Ontario 
M3J 2E5 

(416) 630-5751 
TWX: 610-492-1337 

DENMARK 

NATIONAL SEMICONDUCTOR DENMARK 

Nyhavn 69 

1051 Copenhagen 

Telephone: (1) 153110 

Telex: 160 39 

♦Microprocessor System Specialist Available 



FLORIDA 

♦CARIBBEAN REGIONAL OFFICE 
2721 South Bayshore Drive, Suite 121 
Miami, Florida 33133 
(305) 446-8309 
TWX: 810-848-9725 

ILLINOIS 

NATIONAL SEMICONDUCTOR 
WEST-CENTRAL REGIONAL OFFICE 
800 E. Northwest Highway, Suite 203 
Mt. Prospect, Illinois 60056 

(312) 394-8040 
TWX: 910-689-3346 

INDIANA 

NATIONAL SEMICONDUCTOR 
NORTH-CENTRAL REGIONAL OFFICE 
P.O. Box 40073 
Indianapolis, Indiana 46240 
(317) 255-5822 

MARYLAND 

CAPITAL REGIONAL SALES OFFICE 
300 Hospital Drive, No. 232 
Glen Burnie, Maryland 21061 
(301) 760-5220 
TWX: 710-861-0519 

MASSACHUSETTS 

♦NORTH-EAST REGIONAL OFFICE 

#8 WallisCt. 

Lexington, Massachusetts 02173 
Telephone: (617) 861-6090 
TWX: 710-332-0166 

MICHIGAN 

♦REGIONAL OFFICE 
23629 Liberty Street 
Farmington, Michigan 48024 

(313) 477-0400 
TWX: 810-242-2902 



ENGLAND 

NATIONAL SEMICONDUCTOR (UK) LTD. 
The Precinct 

Broxbourne, Hertfordshire 

EN 107 HY 

England 

Telephone: Hoddesdon 69571 
Telex: 267-204 

FRANCE 

NATIONAL SEMICONDUCTOR FRANCE 
EXPANSION 10000 
28 rue de la Redoute 
92-260 Fountenay Aux Roses 
Telephone: 660.81.40 
Telex: NSF 25956F+ 

GERMANY 

NATIONAL SEMICONDUCTOR GmbH 
8000 Munchen 81 
Cosimastr. 4/1 
Telephone: 089/915027 
Telex: 05-22772 

HONG KONG 

NS ELECTRONICS (HONG KONG) Ltd. 

11th Floor 

4 Hing Yip Street 

Kwun Tong 

Kowloon, Hong Kong 

Telephone: 3-411241-8 

Telex: 73866 NSE HK HX 

Cable: NATSEMI 



MINNESOTA 

DISTRICT SALES OFFICE 

8053 Bloomington Freeway, Suite 101 

Minneapolis, Minnesota 55420 

(612) 888-3060 

Telex: 290-766 

NEW JERSEY 

DISTRICT SALES OFFICE 

140 Sylvan Avenue 

Englewood Cliffs, New Jersey 07632 

(201) 461-5959 

TWX: 710-991-9734 

AREA OFFICE 

14 Commerce Drive 

Cranford, New Jersey 07016 

(201) 272-3344 

TWX: 710-996-5803 

NEW YORK 

CAN-AM REGIONAL SALES OFFICE 
104 Pickard Drive 
Syracuse, New York 13211 
(315) 455-5858 

OHIO 

DISTRICT SALES OFFICE 
Financial South Building 
5335 Far Hills, Suite 214 
Dayton, Ohio 45429 
(513) 434-0097 
TWX: 810-459-1615 

TEXAS 

♦SOUTH-CENTRAL REGIONAL OFFICE 
5925 Forest Lane, Suite 205 
Dallas, Texas 75230 
(214) 233-6801 
TWX: 910-860-5091 

WASHINGTON 

DISTRICT OFFICE 
300 120th Avenue N.E. 
Building 2, Suite 205 
Bellevue, Washington 98005 
(206) 454-4600 



ITALY 

NATIONAL SEMICONDUCTOR ITALY 
Via Valassina 24 , 
20159 Milano 
Telephone: (02) 688 4617 
Telex: 36-540 

JAPAN 

♦NATIONAL SEMICONDUCTOR JAPAN 
Nakazawa Building 
1-19 Yotsuya, Shinjuku-Ku 160 
Tokyo, Japan 
Telephone: 03-359-4571 
Telex: J 28592 

SWEDEN 

NATIONAL SEMICONDUCTOR SWEDEN 
Sikvagen 17 

13500 Tyreso-Stockholm 
Telephone: 08/7 1204 80 
Telex: 11293 

TAIWAN 

NS ELECTRONICS (HK) LTD. 
TAIWAN LIAISON OFHCE 
#60 Teh Hwei Street 
P.O. Box 68-332 
Taipei Taiwan ROC 
Telephone: 563354 
Cable: NSTW TAIPEI 
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